Synchronization in the non-extensive Kuramoto model 
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■ We study the onset of synchronization in square lattices of hmit cycle oscillators with long-range 

P5 ' coupling by means of numerical simulations of the Kuramoto model. In this regime the critical 

^ , coupling strength depends on the system size and interaction range reflecting the non extensive 

' ' ' behavior of the system, but an adequate scaling removes the dependency and collapses the long- 

1/^ \ range synchronization curves with the one resulting from a system with uniform coupling. 

r— I. 05.45.Xt,05.70.Fh 
^ ■ 

^ \ Coupled limit-cycle oscillators have been extensively used to model the collective behavior of systems in physics , 

a chemistry and biological sciences One of the simplest models to describe those systems is the Kuramoto model 

I for the time evolution for the phases of the oscillators. As shown in Q, it possible to reduce any system of N weakly 
coupled limit cycle oscillators near the Hopf bifurcation, to the set of equations 
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where 9i and uji are the phase and natural frequency of the ith oscillator respectively. The coupling term Tij(9i — 6j) 
is a periodic function that may be obtained from the original equations of motion, regardless of the dimensionality 
of the oscillating field. The case of uniform, or mean field, coupling has been considered with detail by Kuramoto [J] 



^ ' for the simplest form of the coupling term, 

T,,^N-^Ksm{ej-e,), iK>0) (2) 



where the factor appears in order to ensure that the typical strength of the net local field experienced by each 
oscillator is independent of the total number N. Kuramoto considered the possibility of a collective oscillatory motion 
■ when the natural frequencies are picked from a given distribution g{uj), as would happen in a real system. As the 
\^ coupling modifies the frequencies, acting as an external oscillatory force, it is necessary to define asymptotic frequency 
of the ith oscillator after the transient regime as 
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, The existence of a collective behavior with phase and frequency entrainment depends on the coupling strength and on 
I ' the choice of g{uj) |4|j|]. In what follows we consider only the case where g{uj) is Gaussian. Kuramoto characterized 
' O \ the degree of synchronization in the asymptotic regime by means of two order parameter like quantities: R for the 
frequencies, and a for the phases. Frequency entrainment may be quantified simply by analyzing the formation of 
frequency clusters and defining 



R= lini ^, (4) 



^ r where Ns is the size of the largest cluster of oscillators with a common asymptotic frequency uj. The complex quantity 



1 ^ 

cr = — ^exp(i6lj) (5) 



is appropriate for the phase ordering as it measures the concentration of phases on a certain value. R and a behave as 
equilibrium order parameters in the sense that, for K below some critical value Kc, oscillations are incoherent leading 
to R and a—Q. For strong enough coupling, or K > Kc, collective oscillation takes place and R,a > 0. 

For the case of uniform coupling, phase and frequency entrainment appear simultaneously at the critical value 
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Kc = (6) 

where ojq is the average value of the distribution of natural frequencies 0| . 

The synchronization in lattices of oscillators with nearest neighbor coupling has also been studied with equilibrium 
statistical mechanics tools. The existence of a critical dimension for synchronization was questioned in as macro- 
scopic synchronization would require R to be finite in the limit N ^ oo. This issue was clarified analytically in |^ 
and numerically in and The behavior of the synchronization curve R x K for different systems sizes in two 
dimensional lattices clearly showed that there is a limiting curve corresponding to the thermodynamical limit 
and, although far from equilibrium, the system is extensive. As opposed to the mean field case, there is no need for 
correction factors in the coupling term. 

Here we focus on systems with long range distance dependent coupling of the form 

T,j = ^smie,-9,), K>0, (7) 

where r^j is the distance between oscillators i and j. The limit a = corresponds to uniform coupling, and a = oo 
to nearest neighbor coupling. From the results obtained in those limits, one can expect that there is a value of a 
below which the system is non-extensive due to the long range interactions. In Hamiltonian systems the limit of non- 
extensivity can be found by requiring the upper value of the energy to be finite as — > cx). Limit cycle oscillators are 
open systems and energy is a meaningless quantity so, instead, we establish the finiteness condition for (l/N) 9i 
in the asymptotic regime. This condition may be established by requiring that the upper value of the coupling term 
divided by N must be finite, as in the mean field coupling case. With this, we end up with a condition identical 
to the one found for extensivity of the Ising model with distance dependent coupling ||9|]l0|, for a < d the system is 
non-extensive. For any value of d and a, it is possible to incorporate the dependence on a and N of the coupling 
term into the value of K with the definition of a scaled coupling strength 

K*{a,N)^KN\ (8) 

where 

N*{a,N)= drr'^-\-'' = — (9) 
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For a > d provides a correction for finite size effects, and for a < d it removes the dependency of the critical 
coupling strength on a and N . For the mean-field coupling {a = 0) and N* = N , for a = d N* behaves like In TV. 

In order to numerically study the non-extensive region, we used two dimensional square lattices of coupled oscillators 
with natural frequencies sorted from the Gaussian distribution g{uj) = (l/\/27r) exp(— a;^/2). Equation (|]) with the 
coupling (0) was then integrated using the Euler method with d< = 0.1. The first 3000 iterations where discarded as a 
transient, and after 6000 iterations past the transient, we calculated the asymptotic frequencies as defined in (^). We 
used the Hoshen-Kopelman method [pi]] to analyze the frequency pattern, two sites were considered to belong to same 
cluster if the frequency difference between them was less than 2i: /T , where T is the time elapsed after the transient. 
R was then associated to the largest cluster, and averaged over 10 initial conditions for the natural frequencies. 

The non-extensivity is evident if we consider the behavior of (i?) as a function of K for a < d. For fixed a, the 
critical value of Kc decreases with increasing N as can be seen in Fig. (|l|.a), as more terms contribute to the coupling. 
The shift in Kc may be corrected with the use of (|8|). In fact, if we plot (i?) as a function of K* , there is no dependence 
on the system size, the curves collapse and a unique value of K* appears. Figure (|l|.b) shows the behavior of {R) as 
a function of K* for fixed a. 

The dependence of Kc on a is evident if we plot (i?) for fixed N and different values of a. Figure (§.a) shows 
the curves for a = 0.25, 0.50, 0.75, and (mean-field coupling without the correction factor Consistently 
with dependence on system size, as a decreases, more terms are included in the coupling, leading to smaller values 
of Kc- Again, fixes the dependence on a providing a perfect collapse even with the mean-field curve as can be 
seen in Figure (g!b). Notice that the value of AT* in Figures (|]) and (||) is in good agreement with (^, our frequency 
distribution leads to K* = 1.59. 

In summary, we have extended the recently proposed equilibrium non-extensive scaling to a system of coupled 
limit-cycle oscillators far from equilibrium, by requiring the finiteness of the long-range coupling term AT/r" . With 
this requirement it was possible to define a universal coupling strength A"*, independent of a and system size, allowing 
the collapse of the curves (A) x K* for different sizes and exponents, including the case of uniform coupling. This 
result allows the study of more general coupling forms, for example to include lattice deformation effects, for any 
range of the coupling term. 
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FIG. 1. Frequency entrainment in square lattices with A'' oscillators interacting through a long range coupling T^/r"'^^. (a) 
The non-extensivity appears in the A'^ dependency of K^, the critical value for synchronization, (b) The corrected coupling 
intensity, K* , absorbs the non-extensive part of the coupling term leading to a synchronization curves independent of the 
system size. 



3 




FIG. 2. Frequency entrainment in square lattices with 1024 oscillators interacting through a long range coupling K/r" . (a) 
The non-extensivity appears in the a dependency of Kc, the critical value for synchronization, (b) The corrected coupling 
intensity, K* , absorbs the non-extensive part of the coupling term leading to synchronization curves independent of a The 
curve for a = correspond to uniform or mean-field coupling. 
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